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A  snace-based ,  ground- looking  radar  must  contend  with  a  propagation 
medium  characterized  bv  properties  that  can  vary  substantially  alone  the 
path  of  propagation.  In  order  to  realistically  predict  the  performance  of 
such  a  radar  one  must  be  cognizant  of  the  properties  of  the  non-homoceneous 
propagation  path.  Since  cany  of  these  properties  are  frequency  dependent 
with  different  properties  exhibiting  different  frequency  dependencies,  it 
nay  be  necessary  to  perform  trade-off  studies  in  order  to  select  the  optimum 
frequency  with  respect  to  the  system  requirements  and  the  propagation 
phenomena. 

In  the  following,  a  survey  is  made  of  the  electromagnetic  propagation 
effects  associates  with  the  varying  properties  of  the  lower  atmosphere  and 
ionosphere  which  are  significant  for  earth-space  propagation  in  the  0. 2-3.0 
GHz  frequency  r?  age.  The  effects  of  weather-related  phenomena  are  included 
in  the  survey. 

2.  THE  LOWER  ATMOSPHERE 

In  this  section  the  propagation  properties  of  the  portion  of  the 
atmosphere  below  an  altitude  of  SO  kn  are  examined.  (Propagation  workers 
often  refer  to  this  lower  atmospheric  region  as  the  troposphere,  although 
in  the  usual  me tero logical  terminology  the  troposphere  is  the  portion  of 
the  atmosphere  extending  from  around  level  to  about  10  km  in  altitude, 
while  the  portion  of  the  atmosphere  lying  between  about  10  km  and  50  km  in 
height  is  called  the  stratosphere  ([l],  op.  953.  1040).)  The  physical 
phenomena  of  the  lower  atmosphere  that  affect  propagation  are  the  gradients 
in  the  refractive  index,  the  scattering  and  absorption  of  electromagnetic 
waves  by  weather-related  precipitates  and  dispersions  (rain,  snow,  clouds, 
fog),  and  the  absorption  and  radiation  of  electromagnetic  energy  by 
atmospheric  gases.  The  consequences  that  these  phenomena  have  for  propaga¬ 
tion  are  considered  in  turn. 

2.1  Refractive  Index  Effects 


The  refractive  index  of  the  atmosphere  is  approximately  l .0003  at  sea 
level  and  decreases  toward  unity  with  increases  in  altitude.  The  gradient 
in  the  refractive  index  can  cause  propagating  rays  to  bend;  the  perceived 
direction  of  a  radar  target  as  determined  from  the  refracted  rays  can  thus 
be  different  from  the  actual  direction,  thereby  introducing  errors.  Doppler 
errors  can  also  be  introduced  if  the  rays  between  the  source  and  the  target; 
propagate  along  curved  rathe:  than  straight  paths.  Finally,  range  errors 
can  be  introduced  both  by  curvature  ol  the  ray  paths  and  the  non-unity 
value  of  the  refractive  index.  Some  quantitive  results  on  these  effects 
will  now  be  presented;  the  angular  error  will  be  considered  first. 

In  {2]  a  graph  is  exhibited  snowing  the  angular  error  due  to  refraction 
when  a  target  in  the  lower  atmosphere  is  observed  from  the  ground.  Employing 
the  data  from  this  graph  one  finds  that  the  refraction  error  at  the  ground , 
&Oq,  for  a  target  at  l 00k  ft  altitude  is  approximately  given  by 
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0.00020  ctn3  for  100%  humidity 

('.) 

0.00015  crr.5  for  0%  humidity 


in  a  standard  atmosphere  (defined  in  ill,  o.  61 )  for  elevation  ancle  5 
greater  than  5°,  where  Aa^  is  in  radians.  In  the  case  of  interest  here  the 
observer  would  be  elevated  and  the  target  near  the  ground,  so  that  this  for¬ 
mula  cannot  be  applied  directly.  However,  by  using  a  formula  for  che 
total  bending  ~  of  a  ray  propagating  through  the  lover  atmosphere  ( [  3  ] , 
p.  24—5) 

T  =  (n,~l  >ctn3 

where  is  the  refractive  index  at  the  ground  (where  the  target  is  assumed 
to  be  located),  and  by  noting  that  (see  Fig.  1} 


T  =  Aa q  +  Aag 

where  ic0  is  the  angular  error  as  seen  by  the  elevated  observer,  one 
finds  that 

Aa0  =  T  -  Acte  , 


or 


Aa0(8  ) 


(nG  -  1 . 00020 }ctn3  for  100%  humidity 


(n q  -  1.0001 3  )ctn9  for  0%  humidity 


(2) 


If  one  employs  the  values,  obtained  from  formulas  in  {2]  at  zero  altitude,  of 
rvj  =  1.000333  and  1.000262  for  100%  and  0%  relative  humidity,  respectively, 
then  He.  (2)  becomes 


0.000133  ctr.3  for  100%  humidity 


(3) 


0.000112  ctnS  for  0%  humidity  . 


At  0  =  5°,  the  minimum  elevation  for  which  Eq.  (3)  is  valid,  one  obtains  a 
value  of  1.6  nilliradians  for  Aag  for  100%  humidity.  For  higher  eleva¬ 
tion  angles  or  for  lower  humidity  the  angular  error  will  be  less. 


Sq.  (3)  applies  when  the  observer  is  at  an  altitude  of  100k  ft. 
Since  the  interest  here  is  on  space  applications  in  which  the  observer 
would  be  many  times  this  altitude,  a  formula  more  appropriate  for  current 
needs  is  required.  If  Aa  represents  the  angular  error  for  observations  of 
the  ground,  then  it  is  shown  in  Appendix  A  that 


2 


AaCn,S)=A<s0(5} 


[sin20+2{ho/ro)  +  Ch0/r0)2]1/2-  sinQ 


[sin-S-5-2{h/rQ)  ~  (h/r  )2]V2_  sine 


where  h0  =  100k  ft,  where  h  is  the  altitude  of  zhe  observer  Cn>n0),  and  where 
r0  is  the  radius  of  the  earth.  The  angular  error  Acs  is  seen  to  decrease  as 
the  observer's  altitude  increases  with  5  reaaining  constant.  In  terns 
of  the  straight  line  distance  D { h , 5 )  between  the  target  and  the 
observer  one  obtains 

3*.  Ho  •  ® ) 

Aa(h,9)  =  D(h, 6 )  Aa0(S).  (5) 


It  can  be  seen  from  this  last  formula  that  the  error  in  the  calculated  ocsi- 
tion  of  the  target  on  the  ground  due  to  refraction,  proportional  to 
Ac{h,9)x  D(h,S),  will  retrain  constant  as  the  altitude  increases 
for  a  constant  elevation  angle. 

Eqs.  (4)  and  (5)  are  derived  on  the  assumption  that  the  index  of 
refraction  is  unity  above  100k  ft;  the  ionosohere  is  not  taken  into  account. 
The  equations  allow  one  to  gauge  the  magnitude  of  the  angular  error  caused  by 
refraction  in  the  lower  atmosphere,  even  when  ionospheric  refraction  is  pre- 


Results  for  the  effect  of  refraction  in  the  lower  atmosphere  on  the 
doppler  frequency  will  now  be  given.  The  problea  was  dealt  with  in  E  2  j  for 
the  case  in  which  the  observer  was  fixed  on  the  ground  and  the  moving  tar¬ 
get  was  elevated.  The  results  given  by  Eq.  (39)  in  [2]  can  be  applied  to 
the  current  problem  by  replacing  Aoy  with  Aa G  so  that  cne  has 

2fv 

if d  =  -  c  Aa<j(3)cos?  sin£  . 

Here  f  is  the  r.f.  frequency  of  the  radar,  c  is  the  speed  of  light,  Af^  is 
the  refraction-induced  doppler  error,  and  V  is  interpreted  as  the  total 
relative  velocity  of  the  target  with  respect  to  the  observer  (the  observer, 
a  satellite,  will  be  in  notion).  The  quantity  $  is  the  angle  between  the 
(relative)  velocity  vector  of  the  target  and  the  straight  line  drawn  from  the 
observer  to  the  target;  the  cos?  factor  was  mistakenly  dropped  in  12]. 
Since  the  doppler  frequency  fg  is  given  by 

2fV 

fd  =  ~  c  cos^  ' 
one  sees  that 


(6) 


^c  =  AaG(S)sin$  . 

-c 

Substituting  the  expression  for  Ao^  given  by  Eg.  (1)  into  Eg.  (6),  one  obtains 
f  0-00020ctn5  sin*>  for  100%  humidity 

(7) 

9.0001 5ctn3  sinp  for  0%  huniaity. 

At  S  =  53,  the  lowest  elevation  angle  for  which  Eg.  (7)  can  be  employed ,  the 
relative  dcppler  error  is  found  to  be  0.23%  near  y  =90°  for  100%  huniaity.  At 
higher  elevation  angles  or  at  lesser  angles  between  the  velocity  vector  and  the 
observer- to- target  line  the  error  will  be  less. 


The  range  errors  introduced  by  refraction  and  propagation  velocity 
changes  in  the  lower  atmosphere  can  also  be  found  in  [2].  A  formula  which 
gives  a  reasonably  cood  fit  to  the  graphical  data  in  the  reference  for  an 
observer  at  an  altitude  of  at  least  10 0'<  ft  is 


Ad(S)  = 


3.3  cscv  for  100%  humidity 


(8) 


7.3  cscd  for  0%  humidity. 


AO  being  the  one-way  range  error  in  feet  when  the  observer  is  at  elevation 
angle  3.  This  formula  will  tend  to  break  down  for  elevation  angles  below 
5°. 


As  an  example  of  the  type  of  errors  observed,  the  Eg.  (3)  gives  a  value 
of  95.2  feet  for  the  one-way  range  error  at  5°  elevation  and  100%  humidity. 

The  effects  described  in  this  section  are  frequency  independent.  In 
the  following  sections  frequency  dependent  affects  will  be  discussed. 


2.2  Heather  Effects 


The  effects  of  certain  weather  phenomena  on  the  propagation  of  elec¬ 
tromagnetic  waves  in  the  frequency  range  of  interest  will  now  be  considered. 
In  particular,  the  effects  of  rain,  hail,  snow,  clouds  and  fog  will  be  examin¬ 
ed. 


The  physical  processes  by  which  these  phenomena  affect  propagating 
waves  are  absorption  and  scattering.  In  general,  waves  propagating  through  an 
atmosphere  in  which  these  phenomena  are  present  will  be  attenuated  by  both 
absorption  and  scattering.  However,  in  the  frequency  range  of  interest  here 
the  wavelength  (a  minimum  of  10  ca)  is  much  greater  that  the  linear  dimensions 
of  any  of  the  associated  particles  (raindrops,  fog  droplets,  etc.)  and  attenu¬ 
ation  due  to  scattering  is  insignificant  in  comparison  to  that  due  to  absorp¬ 
tion  ((41,  pp.  672-673).  !«ever  the  less,  scattering  cannot  be  ignored  since 
the  energy  reflected  back  to  the  observer  may  constitute  a  significant  amount 
of  clutter. 
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2.2.1  Attenuation 


(i)  Rain 


Attenuation  by  rain  will  be  the  first  topic  to  be  considered  here.  A 
f  emu  la  which  has  been  found  useful  in  estimating  this  attenuation  is  [5] 

A  =  aRb  (9) 


6.39x10"^ 


'  '"GHz 


<  2.9 


l  4.21X10-^ 
as  0.351f0-158 


■  2-9<fGHz  <  5-4 
'  fGK2  4  8-5  • 


I 


c: 


1 
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A  is  the  one-way  attenuation  in  d3/km,  R  is  the  rain  rate  in  m/hr,  and  f qj-2 
is  the  frecuency  of  propagation  in  GHz.  The  values  of  a  and  b  furnished  with 
Ec.  (9)  apply  to  rain  at  a  temperature  of  0°C;  the  effect  of  temperature  on 
attenuation  will  be  discussed  shortly.  Table  1  provides  values  of  A  for  a 
range  of  rain  rates  and  frequencies.  The  attenuation  is  seen  to  increase 
approximately  as  the  square  of  the  frequency. 

The  values  in  Table  i  apply  to  rain  at  a  temperature  of  0°C.  Attenua¬ 
tion  will  decrease  with  increasing  temperature,  and  since  changes  of  3  to  l 
can  be  observed  over  a  40°C  temperature  span,  a  correction  formula  for  tem¬ 
perature  is  desirable.  If  the  attenuation  is  principally  due  to  absorption, 
a  reasonable  assumption  in  the  present  frequency  range,  then  the  attenuation 
is  governed  by  the  equation  ([3],  p.  24-22) 

A  =  -8.l86H\~^2h{(c^-1)/(ew*2)}  ,  (10) 

where  A  is  again  the  attenuation  in  c3/km,  M  the  liquid  water  content  of 
air  in  g/m3  (a  plot  of  M  versus  rain  rate  is  given  in  Appendix  3),  A 
the  wavelength  in  cm  of  the  propagating  wave,  and  ~-j  is  the  complex 
relative  dielectric  constant  of  water.  3v  using  the  Bebve  formula  for 
e-rf  which  appears  in  {41,  p.  675,  with  seme  sligut  notational  changes, 
viz., 

CWQ  - 

=  l+i  +  Gw-  (11) 

*C3 

where  ew,  e*-»»  and  AX  are  temperature  dependent  constants,  and  by  taking 
note  of  the  facts  that  ewo»ca»  and  in  the  current  frequency  range 

(which  can  be  seen  free  the  tabulations  in  14],  p.  675),  one  can  show  that 


1 


(12) 


^ ' £wo~^v 


,)AA 


A  { £  j-2) 

CS 3  '  WO 


Ler  Xj  (?,  A^)  denote  the  expression  on  the  left-hand  side  oi 


Sc. 


(12)  with  its  iaolicit  temnerature  dependence. 


i.e. . 


X. (T,A  )  =  Inf  —  1 

i  cm  le  +?! 


1  c-  +~)' 

w 


(13) 


Then  by  substituting  the  tabulated  values  of  AA,  etc.  ((4],  p.  675,  where 
the  decimal  ooint  is  in  error  in  the  30°  and  40°  values  for  AA)  into  the 
right-hand  side  of  Sq.  (12)  one  can  obtain  the  values  given  in  Table  2. 
It  can  be  seen  fret  Eq.  (12)  that  the  quantity  Aa:<;  which  appears  in  the 
table  is  independent  of  A^.  The  tabulated  ratio*  K^(T,  Aa:i)/Ki(0,  A^) 
is  the  temperature  correction  factor,  i.e.,  the  factor  by  which  one  multi¬ 
plies  the  0°C  attenuation  (in  dS)  in  order  to  obtain  the  attenuation  at 
ten  per  a  ture  T.  This  factor  can  also  be  seen  to  be  independent  of  A _ . 

^  *  i  -  -.i 

It  sneura  be  noted  that  the  A— independence  observed  here  would  not  be 
ror  wavelengths  much  less  than  the  minimum  wavelengths  of  current 
concern  (10  cm).  Finally,  it  is  noted  that  seme  of  the  values  shewn  in 
Table  2  differ  by  a  few  percent  from  those  calculated  by  Gunn  and  East 
(16],  Table  1). 


TABLE  2.  Temperature  Dependence  of  Water  Absorption  Factor 
and  Attenuation  Correction  Factor 


T  (°C) 

• 

AcaVT'  *ca> 

1 

Xi(T,  Xaa)/Si(0,  A^) 

i I 

o 

-o.no 

1.0  ! 

10 

-0.0713 

0.65 

18 

-0.0546 

0.50 

20 

-0.0509 

0.46 

30 

-0.0388 

0.35 

40 

-0.0309 

0.28 

Hq.  (9)  and  Table  1  give  the  one-way  rain  attenuation  per  ha.  Bain 
models  have  also  been  developed  for  predicting  the  total  one-way  attenuation 
in  rain  for  various  rain  rates  and  for  various  elevation  angles  of  the 
source.  One  of  these  models  is  that  of  Lin  as  given  in  (7]: 


S 


(14) 


_ 2636 

A_  =  aR2  ^ .  2636sin3]  - 

t  s  -R.--6.2-r  - 

3  4-G 


H£-re  A*,  is  the  total  one-way  loss  in  c3,  R-  is  the  five-minute  average  rain- 
rate  in  nn/h,  G  is  the  elevation  of  the  ground  above  sea  level  in  kn,  and 
3  is  the  elevation  of  the  source.  Quantities  a  and  a  are  the  sane  as 
those  given  with  Sc.  (9). 

Plots  of  the  total  one-way  attenuation  versus  rain  rate  for  several 
frequencies  are  given  in  Fic.  2. 

This  section  on  rain  attenuatich  is  concluded  with  Fig-  3  and 
accompanying  Table  3  {both  taken  iron  [7])  which  show  the  U.S.  and  global 
distribution  of  rain  enployed  in  another  rain  attenuation  model,  the  Crane 
nodel. 

TABLE  3  -  Point  Rain  Distribution  Values  (oa/h)  versus  Per¬ 

cent  of  Year  Rain  Rate  Is  Exceeded  ( fron  i~J) 
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(ii'  Snow 

A  formula  for  the  attenuation  of  electromagnetic  waves  or apace ting  in 
falling  snow  is  given  by  Gunn  and  Sast  ([6],  p.  536).  In  the  frequency  range 
of  interest  here  only  the  second  tern  of  the  attenuation  expression  is  sign¬ 
ificant,  allowing  one  to  obtain 

A  =  0.00224R/Acn  r .5) 

for  the  one-way  attenuation  A  in  c3/kn  at  0°C  where  the  rate  at  which  the 
snow  falls  is  equivalent  to  liquid  water  fall  rate  of  R  sa/h.  (As  a  rule 
of  thunb,  the  liquid  water  rate  is  1/10  the  snow  fall  rate  (Cl),  p.  863).} 
A  tabulation  of  attenuation  for  several  values  of  R  and  frequency  is  given 
in  Table  4.  Attenuation  in  snow  is  proportional  to  frequency  in  the  pre¬ 
sent  frequency  range,  whereas  it  nay  be  recalled  that  rain  attenuation 
tends  to  vary  with  the  square  of  the  frequency. 


TABLE  4  -  One-Wav  Attenuation  in  Snow  (c3/kn) 


- i - 

j 

~GHm  ! 

A  (mn/h  of  Liquid  Water) 

! 

1 

0.25 

1.25 

2.5 

1.5 

10. 

2.5 

i  l 

0.2 

3.7x10-6 

1 .9x10-5 

3.7x10-5 

7.5x10-5 

1.5x10-4 

3.7x10-4 

0.5 

9.3x10-6 

4.7x10-5 

9.3x10-5 

1.9x10-4 

3.7x10-4 

9.3x10-4 

1.0 

1.9x10-5 

9.3x10-5 

1.9x10-4 

3.7x10-4 

7.5x10-4 

0.0019 

2.0 

3.7x10-5 

1 .9x10-4 

3.7x10-4 

7.5x10-4 

0.0015 

0.0037 

,0 

5.6x10-5 

2.8x10-4 

5.6x10-4 

0.0011 

0.0022 

0.0056 

As  with  rain,  attenuation  in  snow  exhibits  a  temperature  dependence. 
The  attenuation  in  snow  tends  to  decrease  with  decreases  in  tenperature, 
however.  The  attenuation  will  be  given  by  Eq.  (10}  if  one  replaces 
with  £t,  the  complex  relative  dielectric  constant  of  ice.  It  is  conven¬ 
ient  to  define  an  absorption  factor  for  ice  analogous  to  the  K;  factor 
defined  for  liquid  water  in  (13).  In  particular  one  defines 


Ki(T) 


ei-i 
In)  £,+2 


l 

i 


06} 


Values  for  are  given  in  {6]  (Table  2,  p.  525).  One  can  employ  these 
values  to  obtain  the  tenperature  correction  factor  K£(T)/K^(o)  for  several 
values  of  T;  the  values  of  the  factor  are  given  in  Table  5. 

The  values  of  froa  which  the  tenperature  correction  factors  of  Table 
5  were  computed  were  actually  determined  from  measurements  in  the  3  cm  to  9 
cm  band.  However,  it  appears  likely  that  the  tabulated  values  of  the 
factors  apply  to  the  entire  wavelength  span  of  interest  here  (10  cm  to  1 50 
cm)  for  the  following  reasons.  The  data  in  [8]  indicate  that  the  Debye 
formula,  which  has  the  fora  given  in  Sc.  (11),  adequately  characterises  the 
dielectric  properties  of  ice  when  il  is  a  temperature  dependent  parameter 
assuming  values  of  from  6.6  kn  to  49  km  in  the  -0.1  °C  to  -21  °C  tenperature 
range.  Then  because  the  AX  parameter  is  much  greater  than  the  maximum 
wavelength  of  interest  (150  cm) ,  the  Debye  formula  snows  that  the  dielectric 
properties  of  ice  will  exhibit  little  variation  across  the  entire  band  of 
interest  as  well  as  at  even  snorter  wavelengths.  The  wavelength  invariance 
of  the  measured  dielectric  properties  of  ice  in  the  3  cm  to  9  ca  band  tends 
to  indicate  that  this  is  indeed  the  cas; . 

(iii)  Hail 

The  attenuation  properties  of  hail  would  tend  to  be  extremely  vari¬ 
able,  assuming  that  the  individual  hailstones  would  generally  have  a  liquid 
water  coating.  This  can  be  inferred  froa  the  finding  that  the  absorption 
in  an  ice  sphere  surrounded  by  a  spherical  liquid  water  shell  can  vary 


frca  a  snail  fraction  of  the  absorption  that  would  be  exhibited  by  a  water 
sohere  of  the  sane  total  nass  to  core  than  twice  that  of  the  water  sphere, 
depending  on  the  thickness  of  the  liquid  water  shell  [6] . 


(iv)  Clouds  and  Foe 


Attenuation  in  fog  and  liquid  water  clouds,  where  the  droplet  size  is 
on  the  order  of  0-01  ca  or  less  ([3],  p.  24-22  and  {4],  p.  677),  is  given  by 
3q.  (10)  with  the  aid  of  definition  (13)  and  Table  2.  The  value  of  M, 

TABLE  5.  Temperature  Dependence  of  Ice  Absorption  Factor 
and  Attenuation  Correction  Factor 


T(°C) 

^(T) 

-  I 

3Ci(T)/Ki(0)  1 

l 

0 

-9.6x10-4 

! 

1.0 

-10 

-3.2x10-4 

0.33 

-20 

-2.2x10-4 

0.23 

the  liquid  water  content,  varies  fraa  l  to  2.5  g/n3  in  clouds,  although 
values  as  high  as  4  g/n3  have  been  observed  ([3],  p.  24-22);  in  fog  M  is 
typically  sruch  less  than  l  g/a3  although  values  as  high  as  2.3  g/a3  say 
be  found  in  heavy  sea  fogs  {[4],  op.  677,  683) . 

Eq.  (10)  can  be  enployed  to  obtain  the  attenuation  in  ice  clouds  by 
replacing  £w  with  tj  and  by  using  definition  (16)  and  Table  5.  The 
values  of  attenuation  obtained  apply  to  ice  crystals  uncoated  by  liquid 
water;  a  water  coating  will  cause  the  type  of  variability  in  attenuation 
noted  for  hail.  The  aass  of  ice  per  unit  volume  of  air  in  ice  clouds 
tends  to  be  substantially  less  than  the  mass  of  water  found  in  liquid 
water  clouds,  often  being  less  than  0.1  g/a3  and  rarely  core  than  0.5 
g/n3  ([6],  p.  24-22). 

Plots  of  attenuation  in  fog  and  clouds  at  0°C  are  displayed  in  Fig. 
4.  Attenuation  at  other  temperatures  can  be  obtained  by  multiplying  the 
attenuation  obtained  from  this  figure  by  the  value  in  the  third  column  of 
Table  2  (for  liquid  water)  or  Table  5  (for  ice)  corresponding  to  the 
desired  tempera  tare . 

2.2.2  Scattering 

All  of  the  weather  phenomena  considered  in  the  section  on  attenuation 
will  reflect  electromagnetic  energy  back  towards  the  source.  This  energy 
nav  constitute  a  significant  aaount  of  clutter. 

For  collections  of  particles  whose  diameters  are  such  less  than  the 
wavelength  of  the  scattered  wave,  as  will  be  the  case  here,  the  scattering 
cross-section  n  in  square  asters  per  cubic  meter  of  the  particle-filled 
aediua  is  (13],  p-  24-22) 
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K  *^m  CV^*- 


n  =  io"10-5x  ^IkI^ 


s  -1 

K  =  e-r2  , 


£  being  the  conplex  relative  dielectric  constant  of  the  particle  scatterers, 
and  where  Z  is  the  reflectivity  factor  in  units  of  niilineters  to  the  sixth 
power  per  cubic  cater,  values  for  which  will  be  given  shortly.  The  factor 
10-10  which  appears  in  (17)  bat  not  in  [3]  arises  fron  the  units  selected 
here  for  n  and  the  wavelength. 

The  factor  |x] 2  which  appears  in  Ea.  (17)  is  relatively  insensitive 
to  both  teaperature  and  frequency  for  water  and  ice,  as  a  study  of  Tables  1 
and  2  in  {63  reveals-  One  can  sinply  enploy  the  values 


UP  = 


0.93  for  iicuic  water 


0.21  for  ice 


incv.  endently  cf  teaperature  and  frequency. 

The  reflectivity  factors  for  the  weather  phencnena  cf  interest  here 
are  {{3],  po-  24-31,  24-32) 

200  R1*6  for  rain 

Z  =  1000  R1*6  for  snow  (19) 

4-8x10~2^2  for  f0g  and  liquid  water  cl-oucs 

where  Z  is  in  units  of  (sa*)6/a3r  where  R  is  the  rain  rate  (or  equivalent 
rain  rate  for  snow)  in  ssm/h,  and  where  M  is  the  mass  of  liquid  water  per 
unit  vo lure  cf  air  in  g/m3. 


A  plot  of  the  volumetric  reflectivity  (scattering  cross-section  per 
unit  volune)  for  rain  and  snow  is  shown  in  Fig.  5,  while  Fig-  6  shows  the 
corresponding  plot  for  fog  and  liquid  water  clouds. 

The  scattering  cross-section  of  hail  should  be  strongly  dependent  on 
the  relative  thickness  of  the  water  coating  the  hailstone  is  likely  to 
have.  This  is  due  to  the  fact  that  the  scattering  cross-section  of  a 
spherical  drop  of  liquid  water  is  factor  4.3  greater  than  the  cross-section 
of  the  sane  drop  when  completely  frozen;  a  partially  frozen  drop  of  water 
consisting  of  a  spherical  particle  of  ice  surrounded  by  a  spherical  shell 
of  liquid  water  has  a  scattering  cross-section  which  exceeds  that  of  the 
completely  frozen  drop  by  a  factor  between  unity  and  4.3,  cepe:  iing  on  the 
relative  masses  of  ice  and  liquid  water  [6]. 

The  reflectivity  factor  Z  for  an  ice  cloud  will  be  the  sane  as  that 


for  a  liquid  water  cloud  if  the  relative  distribution  of  particle  diameters 
is  the  same.  In  this  case  the  reflectivity  of  the  ice  cloud  will  be  about 
0.23  the  reflectivity  of  the  liquid  water  cloud  by  virtue  of  the  ratio  of 
the  values  of  !'<J  -  fo.  ice  and  liquid  water  (Eq.  (18)),  assuming  that  the 
ice  crystals  do  not  have  appreciable  liquid  water  coatings.  Liquid  water 
coatings  would  cause  the  same  type  of  variation  in  cross-section  that  was 
noted  in  connection  with  hail. 


2.3  Gaseous  Effects 


2.3.1  absorption 

The  oxygen  molecule  possesses  a  magnetic  dipole  moment  which  causes 
gaseous  oxygen  to  undergo  an  energy-absorbing  interaction  with  any  time- 
varying  electromagnetic  field  that  may  be  present;  similarly,  the  water 
molecule  possesses  af)  electric  dipole  moment  which  likewise  causes  water 
vapor  to  undergo  an  energy-absorbing  interaction  with  a  time  varying  field 
(14],  pp.  648-664;  {3},  pp.  24-13,  24-14;  [6];  [7]).  As  a  result  of  these 
interactions,  an  electromagnetic  wave  will  suffer  attenuation  as  it  propa¬ 
gates  through  the  atmosphere. 

The  absorption  in  c3/km  at  sea  level  (1013  mb  atmospheric  pressure) 
and  at  68°F  can  be  obtained  irrom  ( [ 3 ] ,  pp.  24-13,  24-14) 
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where  Aq  and  AH  q  are  the  absorption  coefficients  for  oxygen  and  water 
2  2 


vapor,  respectively,  where 
Av-j/c  =  0.018, 

A'J2/c  =  0.050, 

Av3/c  =  Av^/c  =  0.0906, 
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ana  where  pv  Is  the  absolute  humidity  in  grans  (of  water  vapor)  per  cubic 
meter  of  air. 

A  number  of  values  calculated  with  the  aid  of  Scs.  (20)  and  (21) 
are  given  in  Table  6.  Table  7,  adapted  fron  [1],  p.  444,  gives  values  of  p 
in  saturated  air  (1G0%  of  relative  hunidity)  at  1013  no  atmospheric  pressure 
for  a  range  of  temperatures ;  a  plot  of  these  values  is  shown  in  Fig.  7.  Hie 
attenuation  at  sea  level  due  to  water  vapor  is  given  approximately  by  the 
value  of  Au^q/Pv  ta-'r-en  fron  Table  6  tines  the  value  of  pv  taken 


fron  Table  7  (or  Fig.  7)  tines  the  relative  hunidity  divided  by  100%.  In 
order  to  obtain  aore  exact  values  for  AH  0  or  Aq  at  other  temperatures  or 


pressures,  one  nust  enploy  correction  factors  for  the  i’J/c  values  and 
the  coefficients  used  in  Eqs.  (20)  and  (21)  (see  13],  ?p.  24-14  ,  24-15, 
24-16). 


TABLE  6.  One-way  Attenuation  in  O2  and  in  Water  Vapor  (per  gran  of  water 
vapor  per  )  at  Sea  Level  (1013  nb)  and  at  63°  ?. 


H 

A-  (d3/ka) 

°2 

AH_,c/?v  (dB/kn  per  g/n3) 

0.2 

0.00074 

0.00000025 

0.3 

0.0014 

0.00000057 

0.4 

0.0022 

0.0000010 

0.5 

0.0029 

0.0000016 

0.6 

0.0034 

0.C000023 

0.8 

0.0042 

0.0000040 

1.0 

0.0047 

0.0000063 

2.0 

0.0057 

0.000025 

3.0 

_ 

0.0060 

0.000057 

An  inspection  of  Table  6,  or  the  figures  which  appear  in  the 
references,  reveals  that  the  attenuation  due  to  oxygen  absorption  and  the 
attenuation  due  to  water  vapor  absorption  both  increase  with  increases  in 
frequency  over  the  entire  frequency  range  of  interest  here.  Although  the 
amounts  of  the  increases  will  be  different  for  other  altitudes  and  tempera¬ 
tures,  the  trend  will  remain  the  saae.  By  calculating  values  of  “H  0 
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fron  Tables  6  and  7  one  finds  than  nhe  annenuanion  due  no  water  vapor  will 
be  almost  insignificant  in  conparison  no  nhe  02  annenuanion  an  the  low 
end  of  nhe  frequency  range,  even  an  high  ahsolune  humidities,  bun  an  the 
upper  end  of  nhe  frequency  range  nhe  vaner  vapor  antenuation  could  be  a 
subsnannial  fracnicn  cf  the  02  annenuanion  for  high  ahsolune  hunidinies- 
Since  nhe  cne-vav  annenuanion  due  no  bonh  oxygen  and  vaner  vapor  is  nhe 

TABLE  7-  Ahsolune  Hunidnny  in  Sanuraned  Air  an  Sea  Level  (1013  nbl 

{Adapned  fron  nhe  Encyclopedia  of  Amospheric  Sciences  and  Asnrogeology, 
edinec  by  S.S.  Fair bridge -  Copyrighn  (c)  1967  by  Reinholc  Publishing  Corp . 

Ail  righns  reserved.) 


T(  °F) 

p  (g/n3) 

V 

T 

P 

V 

30 

4.44 

65 

15.7 

35 

5.43 

70 

18.5 

40 

6.55 

75 

21  .6 

45 

7.S6 

80 

25.3 

50 

9.40 

85 

29.5 

55 

11.2 

90 

34.2 

60 

13.3 

95 

39.6 

sun  of  Aq  and  Ag  n,  one  would  incur  linnle  error  in  caicu lazing  annenuanion 
2  2 

by  nealecning  AH  q  an  the  lower  frecuencies,  bun  an  nhe  high  freouencnes 
2 

for  high  ahsolune  humidities  the  error  could  be  on  nhe  order  of  nens  of 
percenn. 

In  order  no  obnain  nhe  zonal  one-way  annenuanion  frcn  an  elevaned 
source  no  nhe  ground  for  gaseous  absorpnion,  one  deremines  A0  and  Aj»  g  as 

2  2 

funcnions  of  alninude  and  nhen  innegranes  rheir  sum  along  nhe  path  of  propa¬ 
gation.  For  elevation  angles  between  6 °  and  90°  nhe  zonal  one-way  annenua¬ 
nion  in  a3  due  no  gaseous  annenuanion.  A-.,  will  obey  a  law  of  nhe  fora 

A_  =  a0csc9  (22) 

where  aQ  depends  on  frequency  and  nuaidity;  for  elevanion  angles  benween  0° 
and  6°  one  nusr  employ  a  ray  nrace  sethod  for  which  no  sinple  resulr  ap¬ 
pears  no  exist  [7], 

By  eaploying  innerpolanicn  with  nhe  dara  in  [7] ,  one  finds  rhar  for  a 
ground  nenperature  of  68 °F  (20°C)  and  for  a  relative  huaidiry  of  42%  the 
value  of  a0  for  Sq.  (22)  an  3  GHz  is  0.036.  To  obnain  an  approximate 
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value  of  a0  at  other  frequencies  one  can  enploy  a  graph  given  in  [3]  (p. 
24-20),  which  gives  the  attenuation  vs.  frequency  for  a  low-initial-eleva¬ 
tion  angle,  300  mile  path,  for  two  geographical  areas  at  two  tines  of  the 
year.  On  the  assumption  that  aQ  should  scale  as  the  attenuation  values 
taken  from  this  figure  in  [3],  the  curve  for  Bismarck,  Ji.D.  in  February 


was  employed  to  scale  the  0.026  --alue  of  aQ  at  3  GHz;  the  results  can  be 


seen  in  Table  3.  The  Bismarck  curve  was  employed  to  get  the  scale  factors 
because  the  absolute  humiditv  (as  oooosed  to  the  relative  humid i tv)  would 


be  ercoected  to  be  relatively  low  in  winter,  and  the  0.036  value  of  a_  at 

O 


3  GHz  was  obtained  for  a  value  of  absolute  humidity  (7.o  g/m  at  63°?, 
42%  relative  humidity)  which  is  too  low  to  have  a  significant  effect  on 
the  attenuation.  At  much  higher  humidities  (say  20  g/n^  or  more  at  sea 
level)  the  values  of  a0  might  underestimate  the  loss  by  tens  of  percent 
at  the  high  frequencies.  However,  since  attenuation  at  the  low  end  of  the 
frequency  range  is  fairly  insensitive  to  humidity,  the  values  of  aQ  for 
these  frequencies  should  not  vary  signficantly  over  the  observed  range  of 
terrestrial  absolute  humidities.  As  a  final  note  on  the  aG  scaling,  the 
value  giver,  in  Table  8  for  2  GHz  is  about  5%  less  than  the  calculated 
value  in  [7]. 


lABLE  S.  Attenuation  Scale  Factor  and  Attenuation  Coefficient 

a„  (the  Attenuation  in  d3  at  90°  Elevation ;  for  Gaseous 
Absorotion  for  Low  Humiditv 


i  -GHz 


ie  Factor 

i 

*=>  ! 
| 

I 

0.15 

C.0053 

0.25 

0.0091 

0.35 

0.012 

0.43 

0.016 

0.53 

0.019 

0.59 

0.021 

0.68 

0.024 

0.38 

0.032 

1.0 

0.036 

2.3.2  Antenna  Temperature 

Oxygen  and  water  vapor  radiate  r.f.  energy  which  contributes  to  the 


noise  temperature  of  an  antenna  having  a  bean  that  intersects  the  atnosphere. 
The  noise  tenperature  contributions  of  these  two  atmospheric  conpcnents  de¬ 
pend  directly  on  the  respective  attenuation  constants  at  the  frequency  of 
interest  ({3},  ?.  39-5;  [9];  [10]).  In  general,  oxygen  will  contribute 
core  strongly  to  the  noise  tenoerature  than  will  water  vaoor  in  the  fre¬ 
quency  range  of  interest  here;  this  is  particularly  true  at  the  low  end  of 
the  frecuencv  range. 

rig.  8,  a  composite  of  figures  f rom  [9 j  and  fill,  snows  the  noise 
tenperature  due  to  oxygen  and  water  vaoor  of  a  ground-based,  narrow  bean  an¬ 
tenna  for  several  elevation  angles;  the  tenperature  looking  towards  the 
ground  fron  space  would  be  the  sane  provided  that  the  oean  was  sufficiently 
narrow.  The  absolute  hunicity  for  which  the  curves  in  the  figure  were  calcu¬ 
lated  was  assumed  to  decrease  linearly  free  10  g/nJ  at  sea  level  to  rero  at 
an  altitude  of  5  km. 

The  maximum  and  minimum  gala tic  noise  components  would  be  important 
for  a  space-based,  ground-looking  system  if  the  surface  under  observation 
had  a  high  reflection  coefficient  and  mirrored  the  sky,  as  the  sea  surface 


3.  THE  ICWOSrHSlS 

The  ionosphere  is  the  partially  ionised  region  of  the  atmosphere  that 
begins  at  approximately  50  km  altitude  and  is  generally  defined  as  ending  at 
about  1C00  km  altitude.  The  free  electrons  in  this  region  chance  the  permit- 
tity  fron  its  free  space  value;  the  interaction  between  the  earth's  magnetic 
field  and  the  free  electrons  produces  an  anisotropy  in  the  permittivity  of 
the  region  [i],  It  is  these  changes  in  the  permittivity  which  produce  the 
electromagnetic  wave  propagation  phenomena  characteristic  of  the  ionosphere. 

The  electron  density  in  the  ionosphere  varies  diumallv,  geographi¬ 
cally,  seasonally,  with  sunspot  number  and  with  other  solar  phenomena. 
The  total  electron  content,  i.e.,  the  electron  number  density  integrated 
vertically  from  the  bottom  to  top  surface  of  the  ionosphere,  can  vary  by 
two  orders  of  magnitude  depending  on  the  time  and  location  it  rs  measured; 
predictions  of  these  variations  are  made  by  the  Air  Weather  Service  of 
the  U.S.  Air  Force  [14] .  Because  the  large  temporal  and  geographical 
variations  in  the  electron  density  will  produce  correspondingly  large 
changes  in  the  magnitude  of  the  observed  propagation  effects,  in  order  to 
allow  for  a  conservative  analysis  the  effects  will  generally  be  discussed 
for  typical  daytime  conditions  when  the  ionization  has  peaked.  It  should 
be  bora  in  mind  that  the  ionization  can  increase  by  a  factor  of  TO  above 
the  typical  or  average  level,  however. 

The  discussion  of  the  ionospheric  effects  will  be  restricted  to  cases 
in  which  the  observer  is  located  above  ionosphere,  i.e.,  at  an  altitude  of  at 
least  1000  km.  Space-based  radar  systems  are  expected  to  operate  at  least 
this  high  in  order  to  obtain  extensive  ground  coverage,  and  hence  the  depend¬ 
ence  of  the  propagation  phenomena  on  the  location  of  the  source  within  the 
ionosphere  would  not  be  needed.  Furthermore,  the  restriction  does  allow  a 
considerable  simp licet ion  to  be  made  in  the  discussion. 


o.  ■ 


Refraction  Errors 


The  calculation  of  the  refraction  error  introduced  by  the  ionosphere 
has  been  examined  by  numerous  investigators  (e.g.  [21,  [151,  [161).  Many 
existing  results  are  for  the  refraction  errors  for  ground  observations  of 
elevated  sources,  in  which  cases  sene  manipulation  is  required  in  order  to 
get  the  refraction  errors  for  elevated  observers  looking  at  a  ground  targets - 

Fig.  9,  obtained  frets  [21,  shows  the  refraction  error  for  a  ground 
observer  receiving  a  200  MHz  siqn.al  froa  a  source  located  in  the  ionosphere 
for  three  different  elevation  angles.  The  curves  in  the  figure  were  derived 
by  using  a  layered  model  of  the  ionosphere  for  daytime  ionization  levels. 
If  one  lets  represent  the  refraction  error  (in  radians)  at  the 

ground  when  observing  at  source  located  at  an  altitude  of  1000  km,  one 
finds  that  the  equation 

A  OgO)  =  0. 00033c  Zn9  (23) 

fits  the  data  in  the  figure  to  an  accuracy  of  better  than  15%.  However, 
since  ionospheric  refraction  remains  relatively  constant  between  0°  and  5 3 
elevation  [131  or  possibly  nay  exhibit  a  slight  increase  as  5  increases  in 
this  range  [16],  Sq.  (23)  should  not  be  employed  at  less  than  10°  eleva¬ 
tion.  Letting  ic.  represent  the  ionospheric  refraction  error  for  an 
observer  at  1000  km  altitude  looking  at  a  ground  target,  it  is  shown  in 
Appendix  C  that  the  following  relationship  can  be  obtained: 

—  (1-ni)cos6-Aag(S)sin9 

Aao(S)=  .  (24) 

[«.♦  JV  -  code]1'2 

ro 

Here  n-,  is  the  height  of  the  observer,  IG'K)  km,  rQ  is  the  radius  of  the  earth, 
and  nj  is  the  real  part  of  the  ionospheric  refractive  index  at  altitude  hf. 
At  frequencies  of  interest  here  one  has  (neglecting  losses  and  magnetic 
bias)  [151 

12  2 

aI  5  1  ~  ~2  (25) 

where  m  is  the  angular  frequency  of  the  propagating  wave,  and  where  cau,  is 
the  (angular)  plasma  frequency  defined  by 

■> 

ne(h)e~ 

ta2(h)  =  m~F"  .  (26) 

p  C  O 

The  quantities  ne(h),  e,  and  ae  are,  respectively,  the  number  density  of  the 
electrons  at  altitude  h,  the  charge  of  an  electron,  anc  the  mass  of  an  elec¬ 
tron. 


The  results  in  Appendix  A  can  be  used  to  obtain  the  refraction  error 
at  altitude  h>h-j  if  there  is  negligible  bending  of  the  wave  between  alti¬ 
tudes  h-j  and  h.  Letting  A<x(9,h)  represent  the  refraction  error  at  alti¬ 
tude  h  where  9  is  still  the  elevation  angle  of  the  observer  measured  from 
frea  ground  level,  one  has 
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Fig-  9  -  Daytime  ionospheric  refraction  error  for  a  ground  observer  of  a  200 
KHz  elevated  source  for  three  elevation  angles. 

(Froa  G-H-  Hillman,  Atmospheric  Effects  on  VKF  and  UKF  Propagation,  Proc . 
IRE,  £6,  pp.  1492-1501  (195S),  Copyright  ©  1958  IRE  (no-  IEEE),  used  by 
permission.  ) 
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-here  Eq.  (23) 
distance  D(S,h) 
(27)  as 


for  Aoq  has  been  employed.  In  terns  of  the  straight 
between  the  target  and  the  observer,  one  can  write 


ia(S,n)  = 


(1-nr— 0.C0033)cos6  0(6,^) 


D(6,h) 


.  r‘1  ,.-,1/2 

[(l+— J2  -  cos 2S j 


(2S) 


Sqs-  (27)  and  (28)  are  not  valid  for  5<10°  since  Eq.  (23),  employed  in 
obtaining  these  two  equations,  is  not  accurate  below  10°  elevation. 


In  order  to  obtain  concrete  results  froa  Eq.  (28)  one  oust  know  nj,  or 
equivalently  ne(n)  in  view  of  Eq.  (26).  Employing  the  Chapuan  fcrauia  for  the 
electron  density  given  in  (21  with  a  scale  height  H  of  100  ka  and  values  of 
1  .25:x1012/33  and  300  ka  for  the  naxiaua  electron  density  S  and  the  height  h 
of  the  aaxisun  density,  resoectivelv,  one  obtains  n  (h.  )=6.2ri010/a3.  3=olov- 

6  t  ' 

ing  this  value  of  ne  in  Eq.  (26),  and  using  the  resulting  value  of  in 
Eq.  (25),  one  finds  that  1  -  n-  =  6.3xl9-3  at  200  KHz.  Using  this  value 
in  Eq.  (28),  one  obtains  for  200  MHz  the  result 


ia(5,h) 


At  6=10°,  the  ainisua  elevation  angle  for  which  Eq.  (29)  can  be  enpioyed, 
one  finds  that  the  refraction  error  is  approximately  equal  to  -0.44  aillira- 
cians  at  h=h-j ,  where  the  negative  sign  indicates  that  the  propagating  wave 
will  approach  the  elevated  observer  froa  below  the  direct  path  to  the  target, 
rather  than  froa  above  the  path  as  shown  for  Aag  in  Fig.  1 . 

It  should  be  noted  that  the  value  of  ne(h-j)  calculated  froa  the  data 
in  (2]  and  enpioyed  in  obtaining  Eq.  (29)  is  a  factor  of  2  less  than  the 
aeasured  peak  value  of  ne  given  in  [17],  fro a  which  Fig.  10  has  been  taken. 
(‘The  in  this  figure  is  the  total  nunber  of  electrons  in  a  vertical 
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Fig.  10  -  Daytime  ionospheric  electron  density  and  columnar  content. 
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column  of  unit  cress -sectional  area.)  However,  since  the  bending  of  the 
propagation  path  is  still  relatively  slight  at  altitudes  on  the  order  of  h- , 
the  inaccuracies  in  the  electron  densities  at  these  altitudes  are  not  expec¬ 
ted  to  cause  sore  than  a  20^-25%  change  in  the  calculated  value  of  the 
refraction  error,  which  amounts  to  no  sore  than  0. 1  miliiradians  at  3=10°. 
Since  the  purocse  here  is  to  dete mine  the  cagnitude  of  the  propagation 
effects,  this  uncertainty  in  the  refraction  error  is  acceptable.  A  sore 
precise  design  study  aicht  require  a  acre  exacting  analysis,  however. 

One  can  shew  from  Eq.  (29)  Cor  Sc.  27))  that  the  refraction  error 
magnitude  will  decrease  with  increases  in  the  elevation  angle  for  elevation 
angles  greater  than  10°.  The  equation  cannot  be  employed  at  ancles  below 
10°,  as  stated  earlier.  Sc.  (29)  also  reveals  that  the  angular  refraction 
error  magnitude  for  a  given  3  will  decrease  as  observer  altitude  increases 
(since  D ( 3 , h ) >D ( 8 , h , )  for  h>h.).  However,  it  was  noted  earlier  that  the 
equation  is  valid  only  if  there  is  negligible  bending  of  the  propagation 
path  between  altitudes  "n.  and  h.  Since  the  electron  density,  and  conse¬ 
quently  the  refractive  index,  exhibits  a  gradient  even  at  quite  high  alti¬ 
tudes  (see  Fig.  li  taken  from  [17]),  there  will  be  a  bit  of  additional  bend¬ 
ing  of  the  propagation  path  and  the  refraction  error  will  differ  slightly 
free  the  value  given  by  Eq.  (29). 

At  night  the  refraction  errors  will  be  less  than  those  seen  during  the 
day.  In  [2]  the  refraction  errors  seen  by  a  ground  observer  tracking  an  ele¬ 
vated  target  were  found  to  be  a  factor  of  2.7  less  at  night  in  conparison  to 
the  daytime  errors;  the  refraction  errors  seen  by  an  elevated  observer  track¬ 
ing  a  ground  target  should  show  a  sisdlar  decrease. 

The  numerical  results  given  so  far  are  for  a  200  MHz  wave.  The 
ionospheric  refraction  error  in  the  frequency  range  of  interest  here  is 
proportional  to  the  inverse  square  of  the  frequency  [13].  The  error  is  prob¬ 
ably  of  nc  more  than  marginal  importance  at  200  HKz,  and  would  almost  cer¬ 
tainly  be  inconsequential  at  higher  frequencies. 

The  ionosoheric  refraction  errors  were  calculated  without  taking  into 
account  tropospheric  refraction.  Although  in  an  exact  analysis  the  two  types 
of  refraction  errors  must  be  dealt  with  simultaneously  and  cannot  be  separ¬ 
ately  calculated  and  added,  at  the  very  small  refraction  errors  seen  at 
frequencies  of  interest  here  separate  calculation  and  subsequent  addition 
of  the  individual  errors  is  a  very  good  a poroxi nation  [16]. 

3.2  The  Iris  Effect 


The  real  part  of  the  refractive  index  of  a  lossless  or  slightly  lossy 
plasma  will  be  less  than  unity,  and  as  a  consequence  a  wave  incident  on  the 
ionosphere  will  undergo  total  internal  reflection  ( [18] ,  p.  14S)  if  the  ancle 
angle  of  incidence  (angle  between  the  direction  of  propagation  and  the  iono¬ 
spheric  surface  normal)  is  greater  than  a  critical  value  ic-  As  a  cons¬ 
equence,  waves  whose  directions  of  propagation  lie  outside  of  a  particular 
cone  whose  apex  is  the  elevated  observation  point  and  whose  axis  is  normal 
to  the  ionosphere  will  not  penetrate  the  ionosphere;  the  effect  is  known 
as  the  iris  effect  [13]  and  is  illustrated  in  Fig.  12. 
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Fig.  12  -  The  iris  effect. 


Critical  ancle  6^  is  equal  to  sin  ^ (n^)  *133  where  n^  is  the  minimum 
value  of  the  refractive  index  in  the  ionosphere,  i.e. ,  the  refractive  index 
at  the  height  of  maxinun  electron  density.  The  maximum  electron  density 
appears  in  the  F  layer  of  the  ionosphere  at  an  altitude  of  approximately 
300  kn.  Employing  for  <n.J  the  value  cf  2x1012/n3  taken  from  Vic.  10,  and 
utilizing  the  expression  "tor  the  refractive  index  given  by  Sq.  {25),  one 
can  determine  the  critical  angle  for  any  frequency  of  propagation.  Table  9 
lists  tne  vaiues  of  gc  for  a  number  of  frequencies. 

—  ris  effect  could  affect  an  elevated  observer  survevinc  the 
earth's  surface  by  blocking  the  waves  directed  at  certain  portions  of  the 
surface,  rendering  these  portions  unobservable.  However,  one  can  easiiv 
show  that  a  ray  tangent  to  the  earth's  surface  will  intersect  an  ionospher¬ 
ic  layer  at  the  300  km  altitude  at  an  angle  of  about  72.6°  with  respect  to 
the  normal  of  the  layer.  Since  this  value  is  considerably  less  than  the 
smallest  ?c  value  in  Table  9.  there  will  be  no  blockage  here  of  ravs 
by  the  iris  effect  to  any  portion  of  the  earth's  surface.  The  iris  effect 
thus  has  no  impact  on  an  earth  surveillance  system  operating  at  current 
frequencies  cf  interest. 


TABLE  9.  Critical  Angle  of  Incidence  for  the  Iris  Effect 
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fGHz 

1  * 

0.2 

86.4° 

0.5 

88. 6® 

1.0 

89.3° 

2.0 

89.7® 

3.0  | 

89.8°  I 

3.3  Doppler  Errors 


The  bending  of  the  propagation  path  by  the  ionosphere  will  cause  the 
coppler  frequency  inducec  by  target  motion  to  differ  from  what  would  be  ob¬ 
tained  if  the  propagation  path  were  the  straight  line  between  the  target  and 
source.  The  expression  for  the  relative  coppler  error  is  nearly  identical 
to  that  given  by  Eq.  (6)  for  the  error  introduced  by  tropospheric  refrac¬ 
tion,  the  difference  being  the  replacement  AaG(S)  by  A5^(8).  Thus  the 
relstive  coppler  error  caused  by  ionospheric  refraction  is 
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SSS* 


(30) 


^•*~d  ia„(9)sin^ 
fd 


The  expression  for  Aag(S)  for  200  MHz,  1000  km  altitude,  and  day¬ 
time  conditions  was  given  by  Sc.  (23).  Hap  loving  this  in  He.  (30),  one 
obtains 

— G  _  C.00033ctn6sini>  .  l3,) 

£  — 
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Because  of  the  limitations  on  Eq.  (23),  the  above  equation  is  not  valid 
for  3<l0°.  At  altitudes  above  1000  ka  will  differ  slightly  from  the 
value  given  by  Eq.  (23)  due  to  the  additional  ray  bending  introduced  by  the 
longer  propagation  path  in  the  plasna,  as  noted  at  the  end  of  the  section  on 
ionospheric  refraction.  As  also  mentioned  in  the  section  on  refraction,  the 
expression  for  Acl.  nay  be  slightly  inaccurate  due  to  inaccuracies  in  the 
electron  density  model  enployed  in  (2]  at  the  higher  elevations.  These 
inaccuracies  in  the  value  of  Acf^.  will  lead  to  inaccuracies  of  a  similar 
magnitude  in  the  value  of  the  relative  doppler  error  given  by  He.  (31). 

At  night  Aotg  is  a  factor  of  2.7  less  than  its  daytime  value,  as 
mentioned  earlier.  The  relative  doppler  error  at  night  will  then  be  de¬ 
creased  by  the  sane  factor. 

The  previously  mentioned  proportionality  between  the  refraction  error 
and  the  inverse  square  of  the  frequency  implies  that  the  relative  doppler 
error  will  also  decrease  as  the  '  verse  square  of  the  frequency.  Thus  the 
error  at  400  MHz  would  be  one-quarter  of  the  value  given  by  Eq.  ( 31 ) . 

The  total  relative  doppler  error  in  a  system  is  proportional  to  the 
total  refraction  error  introduced  by  ionospheric  and  tropospheric  effects. 
As  mentioned  in  the  section  on  ionospheric  refraction,  the  refraction 
error  for  situations  of  interest  can  be  closely  approximated  by  adding  the 
individual  refraction  errors  due  to  the  ionosphere  and  the  troposphere. 
The  total  doppler  error  can  therefore  be  closely  approximated  by  adding 
the  individual  ionospheric  and  tropospheric  doppler  errors. 

3.4  Range  Errors 


The  ionosphere  is  a  dispersive  medium  which  will  introduce  an  addi¬ 
tional  time  delay  in  a  propagating  signal  over  what  would  be  expected  for 
a  free  space  propagation  path.  This  additional  tine  delay  is  calculated 
from  the  group  velocity  of  electromagnetic  waves  propagating  in  the  medium 
12),  (13). 

Plots  of  range  error  versus  altitude  for  different  frequencies  and 
elevation  angles  can  be  found  in  12).  Employing  the  plot  corresponding  to  a 
decaying  electron  density  above  the  F-layer  distribution  and  a  200  MHz 
frequency,  one  finds  that  the  day  range  error  AD(6)  caused  by  the  iono¬ 
sphere  is  given  quite  accurately  at  1000  km  altitude  by 


ADO)  =  310csc8\ 


wnere 


9  =  cos 
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fcosB - ^ 

U+h,/r  i  , 


(33) 


and  where  h?  =  250  ka;  AD  is  in  units  of  ceters  here.  Pa  rase  ter  6 '  is  the 
angle  a  ray  eaanating  from  the  ground  at  elevation  angle  5  would  sake  with 
the  tangent  plane  at  the  point  the  ray  would  intersect  an  ionospheric  layer 
at  altitude  h->.  Sc.  (32)  is  accurate  for  0>1O°,  and  probably  gives 
cood  results  down  to  8=0°. 


The  plot  froa  which  Sc.  (32)  was  derived  shows  that  the  range  errors 
are  nearly  independent  of  altitude  at  the  700  kn  altitude  and  above,  so  that 
Sc-  (32)  can  be  used  for  altitudes  in  excess  of  1000  ka. 


The  range  error  introduced  by  the  ionosphere  is  proportional  to  the 
inverse  square  of  the  frequency  for  the  frequency  range  of  interest  here, 
just  as  in  the  case  of  the  refraction  error.  At  200  MHz,  Sc.  (32)  gives  a 
value  of  294a;  at  400  MHz  the  error  would  be  cne-fourth  this  value. 


Froa  the  data  in  {2]  one  would  expect  the  night-tine  ?ange  errors  to 
be  less  than  the  daytine  values  by  a  factor  of  about  3.5,  close  to  the  2.7 
factor  found  for  the  refraction  errors. 


3.5  Attenuation 


Below  a  certain  critical  frequency  the  ionosphere  is  essentially 

opaque  to  electronagnetic  waves.  This  frequency  is  given  by  Eq.  (26)  at 

the  altitude  h  of  oaxiaua  electron  density.  Employing  for  the  aaxiaua 

value  of  n  the  2xl  O1^  fc?  fioure  used  in  the  section  on  the  iris  effect, 
e 

one  obtains  12.7  MHz  for  the  critical  frequency,  a  typical  value  for  this 
parameter  and  well  below  the  ainiaun  frequency  of  interest  here. 


Waves  with  frequencies  above  the  critical  frequency  can  propagate  in 
the  ionosphere,  but  they  will  suffer  attenuation  due  to  collisions  between 
the  piasna  electrons  and  the  other  particle  species  comprising  the  iono¬ 
sphere.  This  attenuation  is  frequency  dependent,  being  inversely  propor¬ 
tional  to  the  square  of  the  frequency  in  the  range  of  frequencies  being 
considered  here. 


Fig.  13,  based  on  the  data  in  [2],  shows  the  one-way  daytise 
ionospheric  absorption  loss.  Although  the  data  in  [2]  was  for  a  1000  ka 
altitude,  the  facts  that  the  coluanar  electron  content  is  increasing  very 
slowly  at  this  altitude  and  above  (see  Fig.  10)  and  that  the  electron 
collision  frequency  is  decreasing  with  altitude  ([19],  Table  7.8),  together 
with  knowledge  that  the  attenuation  is  proportional  to  the  integral  over 
the  propagation  path  of  the  product  of  electron  density  and  electron  colli¬ 
sion  frequency,  inplv  that  attenuation  will  increase  very  slowly  above  the 
1000  ka  level.  The  values  in  the  figure  can  therefore  apply  to  observers 
at  several  thousand  kiloneters  altitude. 
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It  can  be  seen  from  Fig.  13  that  the  two-way  ionospheric  absorption 
loss  will  be  less  than  0.4  dS  at  the  minimum  frequency  of  interest  here, 
200  MHz.  The  data  in  [2]  and  f 1 3 5  indicate  that  at  nicht  the  less  will 
decrease  by  a  factor  of  10  to  20. 

Before  concluding  this  section  on  absorption  loss;  two  types  of 
ionospheric  phenomena  which  cause  unusually  high  absorption  will  be 
mentioned.  The  first  of  these  are  the  auroral  absorption  events  which  occur 
in  conjunction  with  luminous  aurora,  mainly  at  the  higher  geomagnetic 
latitudes  1 1 3 3 -  These  events,  which  can  occur  daily  in  the  auroral  zone 
during  the  equinoctial  nonths  and  which  persist  from  several  minutes  to 
several  hours,  can  cause  increases  in  the  absorption  loss  of  frets  6  to  20 
tines  the  typical  daytime  d3  value.  The  second  type  of  ionospheric  phenomena 
is  the  polar-cap  absorption  event,  which  can  increase  the  absorption  loss 
by  a  factor  of  20  to  60  tines  its  typical  daytime  value,  according  to  data 
in  [13].  The  absorption  increase  extends  from  the  polar  caps  down  to  the 
equatorial  limits  of  the  auroral  zones.  The  polar  cap  event  is  caused  by 
an  influx  of  orotons  produced  by  a  solar  flare,  and  is  most  frequently  seen 
during  the  peak  of  the  sunspot  cycle.  The  event  can  persist  for  several 
days  . 

3.6  Polarization  Effects 

The  presence  of  the  earth's  nacnetic  field  in  the  ionospheric  plasma 
results  in  the  plasma  exhibiting  an  anisotropy  in  its  permittivity  [12], 
{191.  Two  nodes  of  electromagnetic  wave  propagation  are  possible  in  this 
anisotropic  medium,  the  "ordinary"  waves  and  the  "extraordinary"  waves,  and 
in  general  a  plane  wave  propagating  in  this  medium  will  be  a  combination  of 
these  two  modes.  Since  the  two  modes  have  different  propagation  constants 
associated  with  then,  a  wave  propagating  through  the  medium  can  exhibit  a 
rotation  of  its  polarization  plane  and/or  a  change  in  the  type  of  polariza¬ 
tion  (e.g. ,  from  linear  to  elliptical). 

For  frequencies  cf  interest  here,  the  type  of  polarization  will  not 
change  during  propagation  through  the  ionosphere  except  when  the  direction 
of  propagation  is  almost  exactly  transverse  to  the  geomagnetic  field  {13}. 
The  transition  angle  6t  between  the  direction  of  propagation  and  the 
geomagnetic  field  at  which  the  type  of  polarization  starts  to  change  is 
listed  in  Table  10  for  several  frequencies  of  interest. 

Although  the  type  of  polarizations  will  not  change  if  the  angle 
between  the  direction  of  propagation  and  the  geomagnetic  field  is  less  than 
the  transition  angle,  the  plane  of  polarization  can  rotate.  This  rotation 
can  cause  a  loss  cf  power  at  the  reception  antenna.  For  example,  for  a 
linearly  polarized  antenna  the  power  loss  would  vary  as  cos2$  where  $  is 
the  angle  between  the  plane  of  polarization  of  the  wave  and  the  plane  of 
polarization  of  the  antenna. 
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TABLE  10.  Transition  Ancle  at  which  Type  of  Polarization  Begins  to  Chance 
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5  =GHz 
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i  0.2 

39.76° 
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j  o.s 

39.90° 

i  i.o 
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39.95° 

1  3.0 

89.98° 
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The  anount  of  polarization  rotation  depends  on  the  angle  between  the 
direction  of  propagation  and  the  geonagnetic  field;  the  rotation  is  a  aazi- 
zua  when  the  direction  of  oropagation  and  the  direction  of  the  geomagnetic 
field  coincide,  i.e. ,  when  the  oropagation  is  strictly  “longitudinal"  [13]. 
Thus  the  anount  by  which  the  plane  of  polarization  rotates  will  depend  not 
only  on  the  elevation  angle  of  the  observer  but  also  on  his  azimuth. 

In  [2]  a  plot  is  given  of  the  phase  chance  between  the  two  propagation 
codes  for  longitudinal  Propagation  of  a  200  MHz  wave  for  a  two-way  path. 
The  polarization  rotation  is  equal  to  one-half  of  the  phase  chance  between 
the  two  components  [131,  so  the  one-way  polarization  rotation  is  equal  to 
one-fourth  of  the  two-way  phase  change,  rig.  14  shows  the  one-way  polariza¬ 
tion  rotation  calculated  frca  the  phase  change  plot  in  12].  The  loss 
scales  in  the  fig-are  show  the  one-way  loss  and  the  two-way  less  (i.e.,  for 
twice  the  one-way  rotation  ancle)  for  linearly  polarized  antennas  due  to 
rotation  of  the  plane  of  polarization.  Although  the  data  in  [2]  was  given 
for  an  altitude  of  1000  km,  the  low  electron  densities  above  this  altitude, 
and  consequently  their  minimal  effect  on  rotation,  allow  one  to  use  the 
data  at  auch  higher  altitudes. 

It  was  noted  earlier  that  polarization  rotation  is  a  maximum  for 
longitudinal  propagation,  and  will  be  less  when  the  direction  of  propagation 
is  not  aligned  with  the  geomagnetic  field.  Thus  the  plots  in  Fig.  14 
represent  upper  bounds  on  the  polarization  rotation  over  all  possible 
azimuths  and  latitudes. 

It  is  interesting  to  note  that  since  both  the  range  error  and  the 
longitudinal  polarization  rotation  are  proportional  to  the  integral  of  the 
electron  density  along  the  path  of  propagation  they  should  both  exhibit  the 
same  dependence  on  elevation  ancle.  Hencz  the  rotation  should  be  propor¬ 
tional  to  esc  3*  where  9*  is  defined  by  (33).  If  one  employs  the  3=90° 
rotation  values  given  in  Fig.  «4  and  coaoutes  the  rotation  at  5=10°  using 
this  assumed  angular  dependence,  the  values  obtained  differ  frca  those  in 
the  figure  by  u%. 

At  night  the  rotation  is  a  factor  3  5  less  than  its  daytime  value  [2j. 
This  is  the  sane  factor  that  was  observeu  in  the  case  of  the  range  errors. 


TWO-WAY  I, OSS  HON  A  LINEARLY 
POLARIZED  ANTENNA  (<IU)  • 


3.7  Scintillation 


Another  icnosoheric  crcoagaticn  phenomenon  which  may  affect  racar 
s vs tec  performance  is  scintillation,  the  temporal  fluctuation  in  anoiituce 


anc  onase  or  waves  vnrcn  nave  traverses  tne  zoncsnnere. 


The  cnencmencn  is 


attributed  to  diffractive  and  focusinc/def ccusinc  effects  of  eioncated 
irregularities  or  inhomogeneities  in  the  electron  density,  primarily  in  the 
F  lave-  [20] ,  [21].  The  irrecu lari ties  typically  have  dimensions  of  i  km 
in  dir  rtions  transverse  to  the  geomagnetic  field  and  lengths  of  4  km  to 
100  km  along  the  field  [21];  they  exhibit  drift  velocities  of  20—300  m/s 
[211. 

Early  data  on  the  scintillation  of  high-frequency  radio  waves  came 
from  radio  astronomical  observations,  as  noted  in  [22],  but  most  of  the 
recent  data  has  come  from  terrestrial  observations  of  satellite  signals 
[23],  [24].  Many  questions  remain  to  he  answered  about  the  phenomenon  and 
consequently  it  is  still  the  subject  of  much  current  research,  as  a  study 
of  the  papers  in  [25]  will  reveal. 

The  geograohical  areas  where  scintillation  is  of  a  sufficient 
magnitude  to  cause  problems  with  radiowave  systems  are  shewn  in  Fig-  15, 
taken  from  [26].  The  density  of  the  shading  in  the  figure  indicates  the 
depth  of  the  fading  associated  with  scintillation,  the  darkest  shading 
indicating  the  deepest  fading.  The  latitudes  shown  in  the  figure  are  the 
geomagnetic  (also  called  invariant)  latitudes,  which  are  based  on  geomag¬ 
netic  field  lines;  the  relationship  between  these  latitudes  and  the  geogra¬ 
phic  latitudes  are  shown  in  Fig.  16,  taken  from  [21]. 

It  is  seen  from  Fig.  15  that  the  regions  where  scintillation  is 
significant  are  the  polar/auroral  regions  and  the  equatorial  region  in  the 
preaidnignt  and  early  morning  hours.  The  physical  processes  believed  to  be 
the  causes  of  the  ionospheric  irregularities  are  different  in  the  equator 
and  polar/  auroral  regions.  In  the  former  the  irregularities  are  attributed 
to  plasma  instabilities,  while  in  the  latter  particle  precipitation  is 
believed  to  be  the  cause  [21  ] .  There  are  likewise  differences  in  the  nature 
of  the  scintillation  observed  in  these  regions  and  in  the  influence  of 
seasonal  and  other  factors  on  scintillation;  data  on  these  subjects  will 
now  be  presented. 

Amplitude  scintillation  measurements  will  be  presented  first.  Since 
the  data  were  taken  under  a  variety  of  different  conditions  one  cannot 
generally  craw  conclusions  from  comparisons  of  the  various  values.  Instead, 
one  should  use  the  data  to  gauge  the  magnitude  of  the  effects  to  be  expected. 
Near  the  equator,  peak- to- peak  amplitude  fluctuations  of  29  cb  have  been 
observed  in  the  259  MHz  signal  from  the  MARIS  AT  I  satellite  in  Natal.  Brazil 
[27],  while  27  cb  fluctuations  were  observed  in  a  1541  MHz  signal  at 
Ascension  Island  curing  a  time  frame  in  which  fluctuations  of  6.3  cb  and  3 
cb  were  recorded  at  Huancayo,  Peru  and  Natal,  respectively  [28].  Amplitude 
fluctuations  of  9  cb  nave  been  observed  at  4  GHz  in  Hong  3Ccng  [29].  In 
the  northern  areas  where  scintillation  is  significant  peak- to- peak  fluctua¬ 
tions  of  30  db  have  been  recorded  at  Millstone  Hill,  HA,  in  the  150  MHz 
signals  received  from  the  U.S.  Navy  Navigational  system  satellites  [21]. 


Data  on  phase  scintillation  are  not  nearly  as  plentiful  in  the 
literature  as  that  on  anolitude  scintillation,  possibly  because  phase  mea¬ 
surements  are  more  difficult  to  make  than  amplitude  measurements .  Phase 
fluctuations  of  a  few  radians  have  been  observed  at  L-banc  (1575  MHz  and 
1228  MHz)  at  Kwajalein,  Marshall  Islands  129].  Fluctuations  of  a  similar 
magnitude  have  also  been  observed  at  Poker  Flats  Pocket  Pap.ce,  just  north 
of  Fairbanks,  Alaska,  in  the  1239  MHz  signal  of  the  multifrequency  beacon 
of  a  Navy  Navigation  Satellite;  simultaneously,  phase  fluctuations  which 
appear  to  have  magnitudes  in  the  tens  of  radians  were  observed  in  the  133 
KHz  signal  of  the  beacon  [30]. 

Knowledge  of  the  fluctuation  frequencies  that  characterize  scintil¬ 
lation  is  also  required  for  an  understanding  of  the  phenomenon.  These 
frequencies  vary  substantially  frca  the  equatorial  to  the  polar /auroral 
regions.  In  the  polar /'auroral  regions  the  period  of  the  peak-to-peak 
flue  :uations  is  typically  1—3  seconds,  while  in  the  equatorial  region  the 
period  is  a  factor  of  2—10  longer  [23].  The  spectra  of  seme  phase  and 
amplitude  scintillation  observed  at  Millstone  Hill,  MA,  can  be  found  in 
[31].  Amplitude  scintillation  spectra  observed  at  Kuancayo,  Peru  at  254 
MHz  and  at  Narssarssuac,  Greenland,  can  also  be  fcunc  in  [32].  In  this 
last  reference,  statistics  on  the  amplitude  scintillation  are  emoioyed  to 
produce  olots  showing  data  channel  reliability  for  different  message 
lengths  and  signal  threshold  levels. 

Scintillation  exhibits  diurnal  and  seasonal  changes,  and  it  is  also 
influenced  by  geomagnetic  and  solar  conditions,  equatorial  scintillation 
tends  to  be  a  maximum  about  the  time  of  the  equinoxes  and  a  minimum  about 
the  time  of  the  solstices  [21],  [23].  Scintillation  increases  with  increas¬ 
ing  sunspot  number  [28],  [33],  [34].  Scintillation  in  the  polar/auroral 

region  increases  during  magnetic  storms  [21]  {disturbances  of  the  geomag¬ 
netic  field  frequently,  but  not  invariably,  caused  by  particle  influx  from 
solar  flares,  and  typically  lasting  3-72  hrs  [35]}.  In  the  equatorial 
region  magnetic  storms  cause  scintillation  to  decrease  during  periods  of 
high  sunspot  number,  cut  the  effect  is  unclear  during  periods  when  the 
number  is  low  [21],  [23]. 

Scintillation  varies  with  azimuth  and  elevation  angle,  the  phenomenon 
becoming  more  pronounced  at  azimuths  parallel  to  the  local  goemacnetic  field 
and  at  elevation  angles  approaching  0°.  Fig.  17,  taken  free  [21],  shows 
the  angular  correction  factors  for  0_,  the  variance  of  the  log  of  the  sig¬ 
nal  amplitude  and  one  of  the  common  measures  of  amplitude  scintillation. 
The  axial  ratios  in  the  Figure  are  the  ratios  of  the  longitudinal  to  tran¬ 
sverse  dimensions  of  the  ionospheric  irregularities  producing  the  scintil¬ 
lations  . 

A  factor  of  particular  interest  which  affects  scintillation  is  the 
frequency.  The  S^  index,  the  variance  of  the  intensity  fluctuations  anc 
another  common  measure  of  amplitude  scintillation  [20],  is  commonly  taken 
to  exhibit  an  dependence  on  frequency,  although  the  -1.5  exponent 

may  be  somewhat  in  error  for  GHz  frequencies  [21],  [33].  The  ms  phase 

(or  coupler)  fluctuations  are  generally  accorded  an  f-'1  frequency  depen¬ 
dence. 


As  a  final  note  on  scintillation,  it  has  been  observed  that  the  phase 
fluctuations  have  a  Gaussian  distribution  whereas  the  intensity  fluctuations 
are  well  characterized  by  the  Hakagani-a  distribution 


p(:<}=  =— 3 - 
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(34) 


a  =  *4  “  , 

S4  being  the  variance  of  the  intensity  fluctuations  nentionec  above;  aore 
elaborate  distribution  functions  which  exhibit  the  observed  correlation  be¬ 
tween  phase  and  amplitude  scintillation  have  also  been  derived  [36]. 


4.  COSCLUSIOHS 


I 


I 


II 


The  najor  trade-off  to  be  aace  in  selecting  a  frequency  in  the  0. 2-3.0 
GHz  range  for  a  space-base  radar,  insofar  as  propagation  and  scattering  ef¬ 
fects  are  concerned,  is  the  relative  insensitivity  of  the  lower  frequencies 
to  weather-related  phenomena  versus  the  relative  immunity  of  the  higher 
frequencies  to  ionospheric  effects.  Clutter  due  to  precipitation  appears 
to  be  the  cost  significant  weather-related  phenomenon;  a  complete  assess¬ 
ment  of  its  effect  on  a  system  requires  a  knowledge  of  the  antenna  pattern 
and  other  system  parameters,  however.  The  most  important  ionospheric 
effect  is  the  Faraday  rotation,  although  scintillation  nay  also  be  a  prob- 
lem.  At  the  lower  frequencies  the  use  of  circular  polarization,  or  some 
means  of  adaptively  changing  the  polarization  of  the  antenna  upon  recep¬ 
tion,  apcears  mandatory;  moreover,  during  extremely  intense  periods  of 
ionization  such  techniques  might  be  needed  even  at  the  higher  frequencies. 
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APPENDIX  A.  Refraction  Error  Above  the  Troposphere  or  Ionosphere 
in  Terms  of  Error  at  the  Upper  Bounding  Surface  of 
the  Medium 

An  expression  will  be  derived  here  relating  the  refraction  error 
seen  by  an  elevated  observer  at  an  altitude  above  the  top  surface  of  the 
troposnera  or  ionosphere  to  the  refraction  error  seen  by  an  observer  at  the 
top  surface  of  the  troposphere  or  ionosphere.  Referring  to  rig.  A-l,  what 
is  desired  here  is  a  relationship  between  5  and  C'»  the  refraction  errors 
at  altitudes h  and  h*,  respectively,  where  hyh'  and  where  h*  is  the  altitude 
cf  the  top  surface  of  the  troposphere  or  ionosphere.  The  propagation 
medium  below  altitude  h*  is  assumed  to  possess  refractive  index  gradients 
so  that  ray  paths  in  this  medium  will  in  general  be  curved.  Above  altitude 
h  *  the  refractive  index  is  assumed  to  be  constant  so  that  ray  paths  are 
straight  lines  in  this  region. 

It  is  assumed  that  the  relationship  between  5*  and  9*  is  known  a 
priori,  along  with  the  facts  that  t, *  is  of  first  order  smallness  (on  the 
order  of  a  few  mi  Hi  radians  for  cases  of  interest  here)  and  that  first  order 
changes  in  9*  result  is  second  order  changes  in  £*. 

Let  D(h,S)  represent  the  straight-line  distance  between  an  elevated 
observer  at  altitude  h  and  a  target  on  the  surface  of  the  earth,  where  3  is 
the  angle  between  the  ground  and  the  straight-line  path  to  the  observer, 
measured  at  the  target.  Referring  again  to  Fig.  A-1 ,  one  has 


OT  =  D(h,5),  (A. la) 


0"T  =  D(h*,5*)  .  (A.lb) 

Now  from  the  law  of  cosines  for  triangle  CTO*  one  has 

CO" 2  =  CT2  *  0*T2  -  2CT  O'T  cos(S»+x/2), 
or 

(r  +h* )2  =  r  2  +  D2(h*,3*)  +  2r .D(h' ,3* )sinc*  .  (A.2) 

o  o  o 

Solving  for  D(h*,9*),  one  finds  that 

D(h*,3* )  =  (r  2sin26'+2r  h ’ +h ' 2) 1 /2-r  sin9 *  .  (A.3) 

o  o  o 

It  will  be  useful  to  know  how  D(h’,8*)  changes  with  small  changes 
in  3*.  Taking  the  derivative  of  D(h,8*)  with  respect  to  S’,  one  finds 
from  Eq.  (A.3)  that 


RAY  PATH 


Fig.  A-l  -  Geoaetry  of  the  refraction  problen. 


d 

as7  D(h*,5*)  =  -D(h* , 3' ) 


_ cos  5  • _ 

[sin23‘+2(h*/r  J+fn’/r  )2]1/2* 
o  o 


(A. 4) 


One  notes  that 

d _  1/2 

|dS*  D(h* # 5* ) I <D{n* , 3* }/{2h*/r0] 


or 


l3D(h',3')l  |d8* |  (A.5) 

D(h*,S* )  <_  1/2 

[2h«/r0] 


By  applying  the  law  of  cosines  to  triangle  CTO  (or  nore  simply  by 
replacing  h*  and  9*  with  h  and  6,  respectively,  in  Sq.  (A.3))  one  can  next 
show  that 


DCn,S)  =  (r  2sin23+2r  h+h2)1/2^  sin3 
o  o  u 


(A.6) 


Angles  ?  and  ?*  can  be  related  by  an  application  of  the  law  of 
sines  to  triangle  TO'O.  One  finds  that 


TO’  P(h*,9*) 

sin?  =  sing*  =  D(h,9)  sin?*  .  (A. 7) 

TO 


Since  ?’  is  very  small  one  can  replace  sin?*  with  ?’.  Furthermore,  it: 
is  clear  from  triangle  TO’O  that 


?  =  ?*  -  (S«-9)<?< 


(A.S) 


hence  ?  is  also  very  snail  and  sin?  can  be  replaced  with  ?.  5q.  (A. 7)  can 

thus  be  written  as 


D(h* ,9 ' ) 

5  =  D(h,8)  C’(3*) 


(A. 9} 


where  the  dependence  of  ?*  on  6'  has  been  explicitly  indicated. 

It  would  be  preferable  if  Eq.  (A. 9)  could  be  written  in  terns 
of  9  alone,  rather  than  in  terns  of  both  9  and  9*.  Let  6?  represent  the 
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cassis 


s*LJ5f-&r  •^MS^vrSv-V* r*« 


chance  in  the  calculated  value  of  5  if  S’  is  replaced  with  5  in  S a  {.-..9 ) . 
One  can  then  write 


D{h,,et)  SDCn'.e*) 

D(h,5)  5C*0’)  *  o(h,3)  5(5’), 


,ag* (S' ) 

L  ^(S') 


OD(h* , 5 * ) }  D(h* , 9* ) 
0(h,,5*)  J  D(h,S • 


-iiliill  oDthJ^S^) 

L  5'(S*)  +  D(h\-9*)  ] 


{A. 10) 


provided  that  6 ’-8  is  of  first  order  saallness.  But  froa  Bq.  {A. 8)  it 


can  be  seen  that 


S'-S  =  C*-'<5* 


hence  the  snallness  of  S*-S  is  established. 


Sow  it  was  stated  earlier  that  5'  will  exhibit  a  second  order  change 
if  S'  chances  by  first  order.  The  quantity  o?*/^’  which  appears  in  3c. 
(A. 10)  will  therefore  be  of  first  order  s tallness.  Next,  froa  inequality 
(A. 5)  it  can  be  seen  that  the  oD/D  quantity  which  appears  in  inequality 
(A.  10)  will  be  on  the  order ^  of  ( 9 • -6 ) / ( 2h  * /rQ ) 1  /"2 .  For  the  troposphere 
h’=20  ailes  and  {rQ/2h*  > '^^lO;  for  the  ionosphere  (r0/2h* )  ' 
has  a  substantially  saaller  value.  Since  5  '-5  is  cn  the  order  of  a 
few  aiiii radians,  the  quantity  (5-8*  )/(2h/rQ)1/2  will  cn  the  order  of 
a  few  hundredths  at  cost.  Hence  froa  Bq.  (A.  10)  one  sees  that  o£/5  will 
be  on  the  order  of  a  few  hundredths  at  cost,  i.e. ,  the  error  incurred  by 
replacing  5*  with  8  in  Bq.  (A. 9)  aaounts  to  a  few  percent  at  cost.  One 
can  thus  write 


0{h*,5) 

orn,S)  5'(9) 


(A. 11 ) 


Xore  explicitly,  by  employing  Bcs.  (A. 3)  and  (A. 6)  one  obtains 


■>2  2  1/2 
. (r0“sin~3+2rQh*+h*  )  1  -r_sinc  . 


r  _  I — — - - - - - 1  r«(3j# 

1  {r02sin~3-r2r0h-*-h^) '  /2-rQsin9 


(A. 12) 


APPE5JDIX  3.  Liquid  Water  Content  of  Rain-Laden  Air 

The  absorption  by  water  droplets  of  an  electromagnetic  wave  propaga¬ 
ting  in  air  depends  upon  the  total  volume  of  water  per  unit  volume  of  air 
rather  than  on  the  shapes  and  size  distribution  of  the  droplets,  provided 
the  droplet  diameters  are  all  nuch  less  than  a  wavelength.  Mow  absorption 
in  rain  is  generally  given  as  a  function  of  the  hourly  rain  rate,  while 
absorption  in  clouds  or  fog  is  given  as  a  function  of  the  mass  (volume)  of 
wafer  oer  unit  volume  of  air.  In  order  to  compare  the  values  for  rain  and 
cloud/fog  absorption  it  would  be  convenient  to  know  the  volume  of  water 
per  unit  volume  of  air  for  various  rainfall  rates.  This  relationship  will 
be  derived  in  the  following. 

To  simplify  natters  let  it  be  assumed  that  the  raindrops  have  discrete 
diameters  d^,  i=l  ,2, . . .  .  Let  R  be  the  rain  rate,  let  n(d^  ;R)  be  the  total 
number  of  raindrops  per  unit  volume  of  air  for  rainfall  rate  R,  and  let  v(d^  ) 
be  the  velocity  at  which  a  raindrop  of  diameter  d£  falls.  Then  one  can  write 

R  =  gi[di3n(dijR)v{di)]  a  (3.1) 


If  W(R)  is  the  volume  of  water  per  unit  volume  of  air  at  rain  rate  R  then  one 
can  also  write 


W{R)  =  |i:[di3n(ci;R)] 
i 


(3.2) 


Let 


n(di;R)  =  H0(R)f(di;R)  (3.3) 

where  NQ(R)  is  the  total  number  of  rain  drops  per  unit  volume  of  air  for  rain 
fall  rate  R  and  where  f(d^;R)  is  the  fraction  of  raindrops  per  unit  volume 
having  diameter  dj_  at  rain  rate  R.  Employing  (3.3)  in  (3.1)  and  (3-2), 
one  finds  that 

H  =  «f  No(H)E[di~f  (di;R)v(di)  j  ,  (3.4) 


1  .3 

**(R)  =  5"-*Ma)E[di  f(di;R)j  . 
i 

Dividing  3q.  (3.5)  by  Eq.  (5.4)  one  readily  obtains  the  result 
Z[  di3f(di;R>j 

W(R)  =  R  5 - : - - 

Eld^fCd^RWtdi)] 


(3.5) 


(3.26) 
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S)j 
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This  is  the  desired  expression  for  the  water  content  of  rain-laden  air  in 
tarms  of  the  rain  rata. 


In  order  to  obtain  concrete  numerical  results  from  Sc.  (3.6)  it  is 
necessary  to  have  a  distribution  f  (c^ ; R)  and  the  velocity  function  v(c{  ).  A 
nv the r  of  distribution  functions  are  available;  the  one  which  will  be  em¬ 
ployed  here  is  the  frequently  used  Lavs  and  Parsons  distribution  as  given 
by  Medhurst*,  from  which  Table  3-1  is  obtained.  The  formula  to  be  used 
for  v(di )  is  obtained  from  Best**  for  the  standard  l.C.A.N.  (International 
Commission  for  Air  Navigation)  atmosphere,  in  particular 


3.32exp(0.0405z){l-ex?[-(d/1 .77) »• *47j  j  t  0.3<c<j5.0 


v(g)  = 


(3.7) 


1.3lex?(0.0290z){l-ex?[-(d/0.316)--754}}#  0.05_<d<0.3 


where  v  is  in  units  of  m/s,  where  d  is  in  units  of  mm,  and  where 
height  above  sea  level  in  km. 


The  results  of  the  calculations  employing  Eqs.  (3.6),  (3.7)  and  the 
Law*s  and  Parsons  distribution  are  given  in  Table  3-2  and  plotted  in  Fig.  3-:. 


*R.G.  Medhurst,  “P-ainfall  attenuation  of  Centimeter  Waves;  Comparison  of 
Theory  and  Measurement" ,  I  ESS  Trans.  Antennas  and  Propagation,  AP-13, 
pp.  550-564,  1965 


**A.C.  Best,  “Empirical  Formulae  for  the  Terminal  Velocity  of  Water 
Props  falling  through  the  Atmosphere",  Quart  .  J.3.  set.  Soc. , 

76,  mm.  302-311,  1950 


APPENDIX  C.  Ionospheric  Refraction  Error  for  an  Elevated  Observer 
in  terns  of  Refraction  Error  at  the  Ground  Target 

Tne  total  bending  if  of  a  ray  is  equal  to  the  sun  of  the  angles  between 
the  ray  and  the  direct  target- to— observer  path  ar  the  ground  and  between  the 
ray  and  the  direct  path  at  the  observer  (see  Fie.  C-1 } : 

y  =  S  -r  $  (C .  *■ ) 

One  also  has  the  following  relationships  between  the  other  parameters  in  Fig. 
C.  1  *: 

e  =  y  -  (a-3),  (C.2> 


a-3  =  (N  -JI'xl O'^ctr.a  =  Anctna  , 
o 


(C.3) 


acosa0= ( a+h ) cosa  ,  (C.4) 


stana _ 

a  =  tana-tana  -  (C.5) 

c 

The  quat  tity  An  is  the  difference  between  the  ground  refractive  index, 
S^xiO-**,  and  the  refractive  index  at  the  observer,  KxiO_t>. 

One  now  employs  Bq.  (C.3)  in  Sq.  (C.2)  in  order  to  obtain 


T  =  e+Anctna  .  (C.6) 

Next,  with  the  aid  of  3a.  (C.4)  one  can  show  that 
acosao 

etna  =  7/2  (C.7) 

[(a-ih)— a2cos2a  ] 
o 


*R.  Weisbrod  and  L.J.  Anderson,  “Single  Methods  for  Coo poring  Tropospheric 
and  Ionospheric  Refraction  Effects  on  Radio  Waves",  ?roc-  IRE,  4700),  pp. 
1770-1777,  1959 


asasaSssEsaaB 


Stoploying  this  in  5a.  (C.6),  one  finds  that 


aAncosa_ 


y  =  £  j- 


o  ■>  o  i/o 
t(a-s-h)_-a“cos-a0]  '  ~ 


(C.8) 


Si 


[f  this  is  in  turn  employed  in  Sq.  (C.  1 ) ,  one  can  obtain  the  result 


9  =  e  + 


a&ncosaQ 

[ (a+h ) 2-a2cos2aQ] 1 /2 


(C.9) 


bI 


Ba.  (C.5)  can  be  solved  for  e  with  the  result 


cfcnu 

=  ®  I  *-  etna 


CC.lO) 


One  then  employs  (C.7)  in  this  equation  in  order  to  obtain 


| 


% 

I 

a 


£  =  5[l- 


asxnao 


[{a+h } 2-a2cos2aQ ] 1 


(C.n> 


When  this  expression  for  e  is  substituded  into  Sq.  (C.9)  one  finds  that 


AncosaQ-5sina0 

[(1+h/a)2-cos2a0]1/2 


CC.12) 


This  is  the  desired  result  for  the  refraction  error  at  the  observer,  9,  in 
terns  of  the  refraction  error  at  the  target,  aQ. 


I 


The  notation  for  the  geoaetric  parameters  used  in  the  nain  body  of 
of  this  report  differ  from  those  enp loved  here.  In  particular,  $  was  repre¬ 
sented  by  Aa^>,  5  by  Aa^,  Oq  by  6,  ana  a  by  rQ.  In  terns  of  these 
parameters  Sq  {C. 12)  assumes  the  form 


AOq  = 


Ancos9-AaGSin6 
[ ( 1 +h/rQ } 2-cos ^9 ] 1 /2 


(C.13) 


